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I . INTRODUCTION 


The problem considered in this paper is to generate auto- 
correlated random variates via computer and study the effects of 
various levels of correlation on goodness of fit problems. This 
is of interest to the Terrestrial Environment Branch, Aerospace 
Environment Division, Space Sciences Laboratory, George C. 

Marshall Space Flight Center, Alabama, and the financial support 
for the project was under NASA contract number NAS8-29286. The 
results will be useful in determining the distribution or esti- 
mating the parameters of populations in which we have correlated 
observations, such as wind speed and temperature. 

The model and simulation techniques are presented in section 
2. Utilizing a Univac 70/46 computer, normal correlated data is 
generated and analyzed in section 3. The Kolmogorov- Smirnov 
statistic and chi-square goodness of fit are used in the analysis. 
Section 4 describes the computer programs and their use with the 
actual programs listed in the appendix. 


2. MODEL AND SIMULATION 


If X and Y are jointly normally distributed/ then it is 
knom that Y|X is normally distributed with mean 
o 

y y + P^- (X - y x ) and variance o 2 y (l - p 2 ). Assume y x - y y • y 

and o x • o y ■ a then 

Y | X - N(y + p (X - y) , c 2 (1 - p 2 ) ) 
or YjX - y + p (X - y) + oVl - p 2 • Z 
where Z ~ N (0,1) . 

Thus it was decided to recursively generate a random sample, 
X l'*** ,X m v ^ a ^ ormu ^ a 

(2.1) X L » y + p (X i _ 1 - u) + al/l - p 2 
Where z i ~ N(0,1) . 

By mathematical induction it can be shown that 

(2.2) E [XjJ ■ y and Var[X^] - a 2 for all i. 

Moreover, with a little patiencs* , the auto-correlation function 
is given by 

(2.3) P xx (t) " 

Thus with (2.1) we are able to generate a sequence of autocor- 
related normal variates. 

After searching the literature it was learned that (2.1) is 
a version of a first order autoregressive process. Mlhram [3] 


2 



defines an autoregressive process X(t) of order m by 

m 

X(t) - y - ZaiiXtt-i) - y ] + e(t) 
x i-1 1 x 

where {e(t)> is a white-noise process of mean zero. Thus we 
could use (2.1) to generate variates of any specified distribution 
with autocorrelation given in (2.3), by choosing the white-noise 
appropriately . 

One other process that might be useful in applications, 
particularly to wind speeds mentioned by Falls II] , is the moving 
average process of order k. This model is of the form 

k 

(2.4) X(t) - y - Ehie(t-i) 

x i-0 

where (e(t)} is a white noise process of mean zero. The advan- 
tage of this model is that the autocorrelation is zero when the 
lag t exceeds k. For example, the second-order moving average 
process 

X(t) - y - e(t) + .75e (t - 1) + .25e(t - 2) 
has autocorrelation function 

C 1 t - 0 

J .577 t - 1 

p xx (T) " < 

xx \ .154 t - 2 

( 0 T > 3 

Thus the correlation goes to zero after each third observation. 

The model given in (2.1) was simulated by first generating k 
uniformly distributed random numbers in (0,1). By the Central 
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k MM 

Limit Theorem it ie known that as k + «, ? • (I X, - k/2)/ Vk/12 

i-1 1 

e 

approaches a standard normal distribution. It is standard pro- 
cedure to choose k ■ 12; thus, for each 12 uniformly distributed 
random numbers we obtained a single realization from a standard 
normal population. By this technique samples of size n - 20, 40, 
100, and 500 were generated. Bach sample, {Z^:i*l, . . . ,n>, was 
transformed via (2.1) to yield our autocorrelated samples 
Cx^ti-l, . . . ,n). Correlations of + .1, + .3, + .5, + .7 and + .9 
were used. 
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3. RESULTS-GOODNESS OF FIT 


To simplify the problem, a mean of zero and standard de- 
viation of one were chosen in all analyses. For each sample 
size, five distinct samples were generated from the N(0,1) 
population. These samples were then checked for fit by the 
Kolmogorov- Smirnov and chi-square statistics before and after 
applying the transformation. An example is presented in Table 
3.1. naturally, for small sample sizes the fit of the initial 
data was not very good in many cases. This was due to the 
fact that the sample mean was not close to zero or the sample 
standard deviation was not close to one. Be that as it may, 
it was observed that after applying the transformation (2.1) 
small correlations did not affect the fit as much as large 
positive correlations. This is pointed out in figures 3.2 
through 3.6. In each instance the probability level of the 
K-S statistic is plotted before and after the transformation. 
The probability before transforming is represented by a dot and 
by a cross after transforming. It is observed that for p - .1 
there ie not much difference in the probabilities; however, for 
p ■ .9 there is a large difference and in many cases the proba- 
bilities approach zero after transforming. Figures 3.7 through 
3.10 illustrate the probabilities for negative correlations. 

It is interesting to note from figure 3.10 that a bad fit was 
corrected by a large negative correlation. The reason for this 
will be given later. 



Standard Deviation = 1.0 Correlation = 0.9 Number * 100 
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The next step in the problem was to generate a sample from 
the N(0,1) population, apply the transformation five time using 
the correlations of .1, .3, .5, .7 and .9 to observe, in ad- 
dition to checking the goodness of fit, how the means and stan- 
dard deviations were affected. An example of the output is 
presented in Table 3.11. 




Mean=0 . 0 

Std.Dev. 

= 1.0 N= 

100 



Normal 

Corr.=0. 1 

Corr.=0. 3 

Corr.=s0.5 

Corr.=0.7 

Corr .=0 . 

Average 

0.0039 

0.0046 

0.0059 

0.0073 

0.0103 

0.0246 

St. Dev. 

0.9030 

0.9023 

0.8934 

0.8695 

0.8263 

0.7392 

K-S Stat. 

0.0766 

0.0933 

0.0794 

0.0688 

0.0900 

0.0990 

Prob . of 

0.6009 

0.3489 

0.5535 

0.7319 

0.3923 

0.2808 

Chi-Sq. 

13.2 

6.2 

17.4 

o 

• 

00 

8.6 

14.2 

DF 

99 

99 

99 

99 

99 

99 


TABLE 3.11 

Theoretically, from (2.2) we would expect the mean and 
standard deviation to remain constant for all levels of corre- 
lation. I’wever, this is not the case since we are sampling 
from n J (0 , 1) population and the sample mean will not neces- 
sarily be zero nor will the sample standard deviation be one. 
From Table 3.11 we see that initially the sample mean is posi- 
tive and after applying increasing levels of correlation the 
sample mean becomes even more positive. Had the initial sample 
maan been negative it would have become more negative as the 
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level of correlation increased. There were five samples of 
size 20, five of size 40, five of size 100 and five of size 
500 generated, and in every case where the sample size was 40 
or greater the same shift in sample means was observed. That 
is, applying a positive correlation tended to shift the sample 
mean away from zero. 

We also observe from Table 3.11 that the standard deviation 
was decreased. Although this did not happen every time, in 
those cases where the standard deviation increased it only in- 
creased slightly. One conjecture that might be made is that 
observations on the tails of the distribution are affected more 
drastically than those close to the mean. More specifically, 
they are shifted towards the mean by the higher correlations 
resulting in a smaller spread. 





Initially, it was believed that negative correlations 
would yield the same results. However, after investigating 
(five samples of size 20, ten of size 40, ten of size 100 and 
five of size 500) it was learned that negative correlations 
tend to shift the sample mean the opposite direction towards 
zero and in a few cases even past zero. An example of this is 
presented in Table 3.12. 



Mean 

=0.0 Std.Dev.=l. 

0 N=100 




Normal 

Corr .=- . 1 

Corr .=- . 3 

Corr. =-.5 

Corr. =-.7 

Corr . 

Average 

0.0680 

0.0616 

0.0496 

0.0381 

0.0262 

0.0097 

St. Dev. 

1.1046 

1.0993 

1.0902 

1.0739 

1.0243 

0.9044 

K-S Stat. 

0.0897 

0.0751 

0.0729 

0.0615 

0.0478 

0.0963 

Prob. of 

0.3963 

0.6260 

0.6629 

0.8434 

0.9763 

0.3122 

Chi-Sq. 

8.8000 

6.2000 

19.2000 

7.2000 

3.8000 

7.2000 

DF 

99 

99 

99 

99 

99 

99 


TABLE 3.12 


As in the case for positive correlations, negative corre- 
lations tend to decrease the standard deviation. However, since 
this did not occur in all cases, no definite conclusion could 
be made. 
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In every case studied a histogram of the data was printed 
and it appeared that in all cases normality was preserved. This 
should be expected since normality is preserved under a linear 
transformation. Thus, the bad fits recorded by the K-S statistic 
was apparently due to the fact that the sample mean and standard 
deviation were not close to zero and one respectively. Since 
negative correlations shift the mean towards zero, this would 
explain why high negative correlation would tend to correct the 
data and give a better fit to N(0,1). 

A possible extension of the problem would be to test the 
goodness of fit by the K-S statistic using X and S for the mean 
and standard deviation and utilize the tables constructed by 
Lilliefors [2] . 

Other extensions of the problem would be to generate non- 
normal data and observe how the goodness of fit is affected by 
correlations or utilize the moving average process of order k 
given in (2.4) with normal and non-normal data. 
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4. COMPUTER PROGRAM USAGE 


The main line program is titled NASAI and utilizes seven 
different subroutines. Subroutine RANDU is a standard random 
number generator. Subroutine KOLMO which uses subroutines 
NDTR and SMIRN computes the Kolmogorov-Smirnov statistic and 
the probability of exceeding the observed value of that statis- 
tic. Subroutine CHISQR computes the chi-square statistic. Sub- 
routine FREQUE calculates frequencies for input for HISTO which 
computes a histogram of the data. 

To use the programs you must supply two data cards: 

CARD 1 

Col 1-10: A large odd integer right justified for 

a seed of the random number generator 
Col 11-14: The size of the sample you wish to 

generate (1 <_ N <_ 9999, right justified) 
Col 15: The number of distinct samples of size 

N you wish to generate d M <_ 9) 

CARD 2 

The data may be transformed via equation 2.1 using 
five distinct correlations. The correlations (positive 
or negative) are placed in columns 5 through 24. Be sure 
to punch the decimal point or use a F4.1 format. 

Col 1 - 4 : blank 

Col 5 - 8 : first correlation 

Col 9 - 12 ; second correlation 



Col 

13 

- 16: 

Col 

17 

- 20: 

Col 

21 

- 24: 

Col 

25 

- 27: 

Col 

28 

- 30: 


third correlation 
fourth correlation 
fifth correlation 

The mean of the population you wish 
to generate your sample from. Punch 
the decimal point or use F3.0. 

The standard deviation of the popu- 
lation you wish to generate your sample 
from. Punch the decimal point or use 


F3.0. 


The output of the program will be a table similar to Table 3.11 
followed by six histograms. The first is a histogram of the 
uncorrelated normal data. Then the next five histograms will 
be of the correlated data with the five different degrees of 
correlation. This procedure will be repeated the number of 
times you specified in column 15 of the first data card. 


21 



i 


Summary 

Several observations can be made from analyzing the data 
generated in this project. First, it should be observed that 
there is no observable distinction between samples of size 100 
and samples of size 500. Thus we may assume that the process 
tends to converge after 100 iterations. Next, the degree of 
correlation built in did adversely affect the distribution. 

The higher the correlation between observations the more the 
data was transformed. This can be observed by investigating 
figures 3.2 through 3.6. The probability of the K-S statistic 
dropped significantly for larger positive correlations. This 
is due to the fact that positive correlations tend to shift the 
sample mean away from the population mean — the larger the cor- 
relation the greater the shift; and, in most instances the stan- 
dard deviation became smaller. This is pointed out by Table 3.11. 
Although no statistical technique was used, it was observed from 
the histograms that normality was preserved, so the apparent lack 
of fit was due to the shift in means and change in the standard 
deviation . 

For negative correlation the changes were not as drastic. 

The shift of the sample mean was towards the population mean 
rather than away as was the case for positive correlations. 

This is probably due to an oscillating effect. This also ex- 
plains why in many cases a bad fit due to a shifted mean was 
corrected to give a better fit. The same change in standard 
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deviation was observed. Also, a study of the histograms for 
the negative correlations indicates that the lack of fit was 
not due to a change in the distribution but rather a change 
in the parameters of the distribution. 
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APPENDIX 
Program Listing* 


P*GBB I***" 


BILM® 
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A F.jMRaH IV ( VL R L43) SOURCE LISTING: 


01/14/ 7b - 1 PAGE 


1 

2 

3 


V 


5 


b 


1 

To (3 

H 

1 

9 


— in 


ii 


12 


T3 


14 


lb 

102 

- 17) 


17 


IP 

18 

IT 


20 


21 

lb 

22 


23 


24 


2b 

99 

26 


27 


2 P 


29 

2 

30 


31 


32 


33 


34 


35 

b 

36 

7 

rr 

3P 

3 

39 


40 


41 


♦ 2 


43 


44 


45 


46 


47 

10 » 

4P 


— w 

8 

bO 



PROGRAM NASA I 
RfcAL^b YFL 

REAL X (b, 1000),KCJ{6),SUMX(6),SUMXSQ(6),STOX(6),STA2(6) 
i ," FR^T 677 Th r l VS V, X uT2 ( 6 ) , 5US JX'( 6 ) , SUSTDx lb") /FTTE ( l4 ) 

DaTA Y,KSTAK, IIMSH,imSPCE/ 'Y', ' •/ 

uko-g 

REAU(S,1) iX,M,M 
r UP-OA 1(110,14,11) 

IK IX.Lt.OSTUP 
R L A |j ( b , 12 )RU, S I ,CHl 
forma r (t>FA.i, 2 ri.o) 

DU 1 02 L = i , 6 
SUSUX ( L ) *0 
SUS r DX ( L ) =0 
CUM I iNllfc 

DOUl. = l,M ' 

WKl IF (6, 18!SI,CHl,N,P(J 

FljRi'iA I ( ’O', • MfcAN=' ,(-■*, 1,5X, 'STD l)fc V* • , Fj . 1 , 5X, 

1 1 N = « , 14/, 1 1X,6< *C0K* t ,F471,2^T1 
DUlbNN= J ,N 
X( 1 , N N ) * 0 . 0 

Dli99K = 1 , UKO ’ 

SUMxSg(K) =0.0 
SDMa (K ) =0. 

SlDA<k)*0 
DlJ 3 J 3 1 , N 

DU2 1*1,1? _ __ 

CALl K AnDU ( 1 X, I Y, YF L ) 

X' 1, J )=X( 1 , J ) + YF L 

X ( 1 , J ) =X ( 1 , J ) “ 6 

DU3k = i , MRU ~ 

If ( J.NE.DGUTlJb 
X(K, 1 ) =X( 1, 1 ) 

GUTU7 ' ' ”” 

X ( K , J ) = S I +RIM K ) * ( X ( K , J - 1 )-Sl l+COl^SDRT ( l-( RO(K )**2 ) )*( X ( 1, J ) ) 
SUMX ( K ) =SUMX ( K ) +X ( K, J ) 

SUMx$u(K)=$UMXSu(K) + X(K,jr**Z 

CUNI IDOL 
DU 1 0 3K = 1 , NRU 
SUMa (K ) =SOMX(K ) /U 

S I Da ( K ) = ( SHMXSU ( K ) - ( N* ( S<DlX ( K ) ) ) ) 

S I Da ( K ) =SUP I ( STltX (K ) / (N-l ) ) 

SUSUX (K ) *SU$UX ( k ) +5UMX (K ) 

SUS IDX(K ) =SUS I DX (K )+STDX(K ) 

L L I. = 1 

CALI KnLMO(X,N,k,STA?,PR?,lLt,ST,CMIVIFin 
CAU CM I SUP ( X,M,K,CMI2,XDF2 ) 

WKl I F ( 0, 8 ) SUMX 

rORMATC ', 'AVFKAGC ,2X,6(F9.5,lX)r 
WKl I L ( 6, 10 ) STUX 



f I ■ ( SnlJ Stt l»j I k: rn ■ 


Okkii 



FjRTRAU IV (VfcK L ^ i ) SOURU LISTING: NASA! PkilbRAM 01/14//b 2 P AGL C 

bl 10 furmatc '.'Sir. • >b< r-9.b> lx n 

b2 WRl ir (b>?b)STA2 

b3 ?b FUP.MAM' ','K-S S I A I •# 6 ( I 9 . •>, 1 X > ) ... 

bA WKI I \. I 6/ A /)PR2 

bb W KURMATC ' > ' PRllB UF ' , 6 ( F9 . b, l X ) ) 

bb WKl IL(6>»W)UI12 _ _ 

b7 87 '" FjJR MAT ( 1 ','CHI-SU ‘ / b ( F 9 . 5/ 1 X ) T 

b8 WKI I L l b / 0 b ) X U F 2 

b9 Bb FURMAK' ' / ' L)t ORE FKF ' > b C F9 . b, 1 X ) / / ) _ _ 

(>n : HI l 1KKK= 1 t b 

ol KRR = Ru ( KK K ) 

t>2 call frlquf i x.kkk^iofrc ) 

b3 CALL HlSTD(RRK/FRt, iA/KSTAR^NSPCti lUAbH) 

b* 1 1 CUN 1 INUL 
bb liUlOAK* 1 .NR[J 

ob Sl'SUX i k ) sVjUsux ( H ) / 1 ~ - 

o7 lviA SOS I MX t K * - SUS I OX ( K ) / M 
00 W|< I I [ ( h, 1 / ) SUSUX.. S'lbTpX 

b9 17 FORMAT! 'O' , 'SUM AVRfc . ' 6 ( F9. b, IX ) W' SlO AVRfc.' ' ' 

7 < > l;b(|9,bslx)> 

a GLJTU1O0 _ 

J> FNO 




O x ^Ji ^ v~n -cM UJ JM 


i 


A F ijP 1 P AM 

1 


1‘1 



| a IW!i»Lf]lkr^« 






v (VfcP hi) bnokct listing: hanl>m subkihjtine oi/14/yb 3 i'agl 

SubKHuTIhE kAN[)U( IX# IY/ yfl ) 

RtAL*a YFc 
I Y = IX*6!>i>39 
' TF f i Y ) 5 , 6,6 
I Y = I Y + 2 i 74 b A<>4 /♦ 1 
YUM Y 

VH = YFL*.Y,6S6o14E-9 

I X= I Y 

RLTUPN 
( NO 


I I 


1 riflfcUffillMIUlllltaA J fc !■>«. 

3 

f i »W I P AM IV (VLI< L4 1) bIJURCi LISTING: KULMU SUIJKIJUT I NE 01/1 »/79 4 PAGE 


1 SUBkUUT 1 NL KULMU(X#N#K#Z#PRUB# I CUO#U,S, I £> ) 

2 t) I MENS I UN X(6#lU0u)#PR0B(6)#Z(6) 

3 ItR = 0 

4 l)U b I=i,N 

5 II (x(K# 1 )-X(K# 1-1 ) ) 1*5,5 

6 1 TEMP = X ( k # 1 ) 

7 ' I M= I - 1 * 

B DO 3 

9 L = I - J 

10 IT ( I EMP-X(K#L> )2#4#4 

11 2 X(K#Lfl )=X(K#L) 

12 3 CUN i iUUt 

13 ' <(K,1)«TCMP ' 

14 Gl)T[J5 

1 54 X(K#L+i)*IEMP 

16 5 CUNiluUt 

17 flHl=N-l 

IB_ X N ~ N 

19 I) N = 0 , 0 

20 Pi» = 0.0 

2! _ I L = 1 _ 

22 6 01J7I *IL,MM1 

23 J=I 

24 It (X(K#J)--X(K#J + 1 ) )v#7#9 

257 CUN 1 INUE 

26 6 J=N 

27 9 IL=J+1 

2« Fl=hS 

29 F S = T LUA1 ( J ) / X N 

30 1 1 ( 1FLHD-2) 10# 13# 1 t 

31 10 II (5 ) lb 11# 12 

32 11 I tR = 1 

33 F,u TO 29 

34 12 ‘ Z(K) = (X(K# J ) -U » / S ' ' 

35 CALL NOIP ( Z#Y#U#K ) 

36 GuTU? 7 

rrn irrsm, 11,14 ~ ' — — 

3B 14 Z (K ) ■ ( X (K# J )-U) /S+l .0 
39 IKZlK) ) I5#15#lt> 

•♦0 Is Y-0.0 

41 GUTU27 

<♦2 16 Y=1 ,-LXP( -Z (K ) ) 

'43 GuTO? 7 

44 W II ( I F C DU-4 ) 1 B # 20# 26 

45 IB II (b) 19, 1 1,19 

4'Tn Y = AfAN((X(K,J»-U)/S)*0.31d3099 + 0.5 

47 GLITu< 1 

4B 20 II ( b~U ) 1 1 # 1 1 # 2 1 

49 M IF(X(k# J )-U)22#22#23 

50 22 Y=0.0 
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A TilPlMAH IV (VLK 14?) bUiJPtL LISTING: KIJLMI) SUUKIIUTINE 

51 GflTlJ?7 

5? 7 3 II (XU, J )-S )75/25/?4 

5* 2* Y-1.0 

~54~ * Gu 11127 

55 ?5 Y=.< a(K/ J >-U)/(S-U) 

56 GljTu2 7 

5TT6 TtR=l 

58 GUTU2V 

59 27 r. 1=AP5 ( Y-I I ) 

6T5 cs*ATib(v-Tsr)' 

61 !>N*AMAXi<liN,EI/fcS) 

62 III 1L-N) 6/8/28 

bT" 2o Z(K)=UN*$0RT(XN) 

64 Call &M1KM 2/PRijK/K ) 

65 ZIK)=Z(K)/SUP KXN) 

66 PRnMM = l .O-PHOrt(k) 

67 2 y i? I TUPN 

I 68 run 



S.H'lPfVi IV < V L f 14 U 3II0KU LISTING: SMlkli SUUkUUT 1 lit 


01/14/76 6 PAGL 


1 suHKiiurniL smikn(x*y#k) 

2 OlMLNSIlili X(6)*Y(b) 

3_ II <X(K)-.27) 1,1,2 

A 1 Y(K ;=<J.U 

3 0UTU9 

hi If < a(K )-l .0) 3,6,6 

7 3 ' Gi = tXf>(-i. 233 / 01 /x(k)** 2 ) 

H G2=G1*G1 

9 G4=G?*G2 _ 

in G o = G 4 * G 4 

11 If <GH-1 .OL-26)4,6,S 

1? 4 Gb=0.0 

13 6 Y(K)s(?.50662 8/A(K) )*Gt*( I.U + g’bVY i .0 + Gfl*G8T) 

14 G(j T U9 

lit 6 I M a ( K ) - 3 . I ) b , 7 . 7 

16 7 Y ( K ) = 1 . u 

17 r.uTuv 

ifln G l=fXiM- 2 . 0 *Xll<)*A(tO) 

19 02=01*01 

20 '.4=02*02 

21 00=04*04 

22 Y(K ) = i .u-,'.U*(Gl-G4 + G«*(Gl-( J H) ’) ----- 

23 9 KLTUPU 

24 HlO 


r 


A FUW1KAN IV (VfcR L<*3) SUURCF LISTING: NUTK SUBRUUTlNE 01/lA/7t> 7 R AGb 


T WftKliuTIhfc NuTr (X* P>f*/K) 

? n illusion x ( 6 ) 

3 X /. » X ( K ) 

A Aa*ABS ( XX ) 

b T - !.()/( l,0+.24lt>AlV*AX) 

b D^O, 39fm*f3*F.XP(-AX*XX/2.U) 

T P = 1 .<)-0*T* ( ((11, 3 3 o 27 <**T - l.B 212 b oT *T+T77FT A 7Tr*T- 

» 10.3!>6!>638)*U0.3l93B15) 

9 IKXX)1,?j2 

irrr p=i.o-p 

11 2 kb TURN 

i? run 



1 


1 


r I jp T « AM [V (VLP l<*3) iUURCL l i S T IN<~, : CHlSok SUkk.lUTlUb Ol/14/7!> 8 RAGb ( 

“T SURkfiUTiNh CHISUkU>n,K,CHlVAL,XI)k ) 

2 OIMtNSIlJM Mb, lU0o>,M9),Cb LL ( 10 ) * X[)F ( 6 ) *CHI VAL (6 ) 

3 I'U 1 1 X= 1, 10 

4 1 CLLL ( IX ) =0.0 

5 CH1VAL(K)-0.0 

6 fl<U«-l.2ttl/ 

7 4( 2)«-0.»»<*ld 

H 013) =-0 . t>244 

9 H ( 4 ) =-0 . 2I>3 J 

10 ' l-Vb) -O.u 

11 >U2 J = (.>» 0 

12 2 3( J ) ■-»( lO-J ) 

lV r/PVAL = HUAl(W)/ 10 . 

14 'Hj3 1 = 1, M 

Lb II ( A(K,l) .1 r ,H( 1 ) )Cbl.L( n = CbU ( 1 ) + l . _ _ _ 

16 ! I ( A ( K . L ) . (, I . >U V ) ) C h 1 L ( 1 0 ) =C F L L ' 10 ) ♦ l” 

L 7 IM 1 i . i= i. , 9 

1R 1K«(k,l) .GI . ) .Atm.A(k,L) .LL.rtlM) >tLLLM)*U Ll.lMl+1, 

”19 3 " cum IijUl 

20 ru6 1 = 1 , lu 

21 s CHlvAL(K)=ClllVAL(K)-»(((CPLL(l ) -b Xk VAL ) »»2 ) / fcXP V Al) 

’ ti. X III (K ) =U- l 

23 Rtluku 

2_4 l : un _ 




I 

1 I A r-JPTHAM IV (VU' L^t 3 ) iUUKtU LISTING: IkFUtib SUHRUUTINE 0 1 / 1 ^ / 7 *» 9 PAG 



6 fk<lk)»hmlk) + i. 

l\ 7 CIJNI1UUL 

2T PlToRn 

& TNI) 







V F.IR1PAM IV (VL« Lk 3) SLIUKU. LIST IMG: HlSIl). SUflRUUT INE 01/14/75 10 PAGt 


1 

2 

3 


<♦ 


s 

l 

6 


7 

4 

n 


9 

It 

10 


11 

s 

1? 


13 

7 

14 


15 


1 6 


17 


ik 

JtO 

19 


20 


21 


22 


23 

h) 

2 A 

'u> 

2 7 

so 

26 


27 


2D 


2D 


jo 


ll 


J? 

7o 

33 


34 

MO 

36 

Mi 

JO 


37 

'no 

3H 


39 


90 


-.1 

Oi) 

<*? 


9 3 

<U 

‘♦A 


’"♦S 



SolmlHjTlUL H l bTl) ( KNU / h Rfc 0 / I N / K / Ni ll H/ MI UUS ) 
DlMLNblU'l JUKI ( l4)/hPEQ< 14)/ 1USM 132) 

DlMttlbTUii B ( 13) 

‘ HI 1 1 = 1, ~ ' “ 

I U S1 1 ( I ) =f’lMUS 

UK III (6/4 ) RKtJ _ __ 

FURHA I ( ' OCPKR= '/ F5 . 1 ) 

Dill 2 1*1/ IN 

JuUl ( I )*FkFu< I ) 

7iYi i c (ft/**) ( anuT i n ) / 11 = 1/ in ) 

F [j k f I A 1 ( ' I RfcfjUL'wCV '/ I4J8) 

‘IK I 11 (6/ 7 ) IUSH _ _ 

FUkllAI ( IX/ I j2Al ) 

F M A X - (J . U 

D Li 201*1/ Ih 

ii ( i k l d ( i ) . l i . r m a x ) on r 1 1 2 o 
FmAa = FRlDI I ) 

Cl)M I 1 I'D) l 

JSC AL * l ~~ " 

II (f MAX.LL.50 .(J)GIjTU40 

JbCAL=(hMAX+49.U>/60.0 

•Ik I If 16/ 30 )K/JS(.AL 

ruRHAM' LACK' /2X/A2/2X/ 'tOUALS ' / I 3/ ' POINTS'/) 
'USD I *1/ III 

.Kill I ( 1 ) = MllTH 

M/,X = rMAA/l LIJA I ( JSLAL) 

'»UB0I = 1/MaX _ 

x* 1 - 1 ) 

D||/U.l = l/ IN 

Ii ( H'LOl J ) /FLUA I (J SLAL ) .1 r .X )0 1 1 1 1 70 
J 1 1 1 1 1 ( j ) = k 

CUN I iNUfc. 

IX = X«f LUA I ( JSCAL ) 

•■.'K I I F (6/ h 1 ) i X/ ( JIM) I U ) / J* i7 IN ) 

F 1 1 K 1 1 A I ( 16 / 5 X / ml 6 x / A ?)) 

dl*>o 1 = 1/ in 

juui i i V* i 

WkI i t ( 6/ 7 ) IUSH 

6< 1 ) =-3.0 

01190517'* 2/ 13 

ill lr.)=f)( I K- 1 ) ♦() . 5 

Wk I IF 16*91) ( III KKK ) / kK K * 1 / 13) 

fckmat:' inTikVal' . iox,n(ri,.UAxn " ~ 

Uhl URN 
l NO 


UMM IUIMI MO. 4 





